Strong Tunneling and Coulomb Blockade in a Single-Electron Transistor 
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We have developed a detailed experimental study of a single-electron transistor in a strong tunneling 
regime. Although weakened by strong charge fluctuations, Coulomb effects were found to persist 
in all samples including one with the effective conductance 8 times higher than the quantum value 
(6.45 kf2) _1 . A good agreement between our experimental data and theoretical results for the strong 
tunneling limit is found. A reliable operation of transistors with conductances 3-4 times larger than 
the quantum value is demonstrated. 



Mesoscopic tunnel junctions between metals represent 
a nontrivial example of a macroscopic quantum system 
with discrete charge states and dissipation |T],|| . Charg- 
ing effects in such systems can be conveniently studied in 
the so-called SET (single electron tunneling) transistors. 
A typical SET transistor consists of a central metallic is- 
land connected to the external leads via two tunnel junc- 
tions with resistances Rl,r and capacitances Cl,r (see 
Fig. 1). In addition to the transport voltage V, the gate 
voltage V g can also be applied to the metallic island via 
the gate capacitance C g . 
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FIG. 1. Schematics of a SET transistor. 

Provided the junction resistances are large Rl.r 3> 
R q = irh/2e 2 ~ 6.45 kf2 tunneling effects are weak and 
can be treated perturbatively || . For a quantitative mea- 
sure of the tunneling strength we define the parameter 
ot% = R q /Ro, where 1/Rq = ^/Rl + ^/Rr- After each 
electron tunneling event the charge of a central island 
changes exactly by e and the energy difference between 
initial and final states of a SET device is typically of or- 
der E c = e 2 /2C, where C = C L + C R + C g . As long 
as tunneling is weak a t < 1 dissipative broadening of 
different charge states T is small (at low T it is roughly 
r ~ a t Ec) and these states are well resolved in energy. 
This ensures nearly perfect quantization of the charge on 
a central island in units of e. As a result at sufficiently 
low T < Ec Coulomb effects dominate the behavior of a 
SET transistor leading to a number of observable effects, 
such as Coulomb blockade of tunneling, modulation of 
the current through a SET transistor by a gate voltage 



V g , Coulomb staircase on the I — V curve etc. ffl. 

The situation changes if the effective resistance Rq be- 
comes of order of R q or smaller, i.e. a t > 1. In this case 
dissipation is large and the excited charge states of the 
system become broadened and overlap. Do strong charge 
fluctuations lead to a complete smearing of Coulomb ef- 
fects in highly conducting mesoscopic tunnel junctions? 

A positive answer on this question was suggested in 
Refs. ||,[|. It was argued that - similarly to Ohmic resis- 
tors (sec e.g. Q) - charging effects should be completely 
destroyed in junctions with a t > 1. According to 
no Coulomb gap can exist for such values of at and a 
crossover to a purely Ohmic behavior was predicted. 

Note that the above conclusion was obtained from 
studies of a certain correlation function which turned 
out to be irrelevant for the problem in question. In con- 
trast, in Ref. a nonperturbative analysis of the junc- 
tion ground state energy has been developed. This anal- 
ysis clearly demonstrates the existence of a nonvanishing 
Coulomb gap E c oc Ec exp(— 2cv t ) in the spectrum of 
the system even for large at 3> 1. Thus strong tunneling 
does not destroy Coulomb effects, it only leads to effec- 
tive renormalization of the junction capacitance. As a 
result the temperature interval relevant for charging ef- 
fects shrinks, but they still remain observable even at 
T > E* c . 

This behavior is qualitatively different from that of an 
Ohmic resistor. The physical reason for this difference 
was also pointed out in (|]. It is due to different sym- 
metries of the allowed charge states: the symmetry is 
continuous in the case of an Ohmic shunt, whereas only 
discrete e-periodic charge states are allowed in the case 
of a normal tunnel junction (see also Ref. 0). The latter 
symmetry remains the same for any SET strength, and 
therefore at low T Coulomb effects survive and can be 
well observed even in highly conducting junctions. The 
results H were supported by subsequent theoretical stud- 
ies (see e.g. (7| for a discussion and further references). 

In spite of all these theoretical developments, an ex- 
perimental investigation of this problem was lacking. Re- 
cently Joyez et al. [pj carried out an experiment aimed 
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to study strong tunneling effects in SET transistors with 
higher junction conductances. Deviations from the stan- 
dard "orthodox" theory Q were clearly detected in four 
different samples [|| . A very recent analysis demon- 
strated that for three samples J8| with smaller a t a re- 
markably good quantitative description of the data || is 
achieved already within the perturbation theory in a t 
if one retains the "cotunneling" terms oc a 2 . Hence, 
the above samples [|| are still well in the perturbative 
weak tunneling regime: the corresponding values of a t 
are large enough to observe deviations from the "ortho- 
dox" theory, but still too small for nonperturbative ef- 
fects to come into play (see below). 

The main goal of the present paper is to develop a 
detailed experimental study of Coulomb effects in meso- 
scopic tunnel junctions in a nonperturbative strong tun- 
neling regime at ^ 1. Beside its fundamental importance 
the problem is also of interest in view of possible applica- 
tions of SET transistors as electrometers. The operating 
frequency range of such devices with at <C 1 is usually re- 
stricted to very low values. One possible way to increase 
such frequencies is to decrease the junction resistance, 
i.e. to fabricate SET transistors with a t > 1. But do 
such SET transistors exhibit charging effects? 

In this paper we investigate if Coulomb effects sur- 
vive in the strong tunneling regime a t > 1, in which 
case discrete charge states are essentially smeared due to 
dissipation. For this purpose we have carried out mea- 
surements of the current-voltage characteristics of several 
SET transistors with various values of the junction resis- 
tance Rl.r ^ Rq in the limit of a zero resistance of the 
external circuit. The results are compared with the ex- 
isting theory plpl . 

Experiment. We have fabricated several SET transis- 
tors with different values of the junction resistance. The 
transistors were made using a standard electron-beam 
lithography with two-layer resist and two-angle shadow 
evaporation of aluminum. Five transistors with junction 
resistances in the range from 2 kfl to 20 kf2 were stud- 
ied. The corresponding values of a* varied between 1.5 
and 8.3. The crossection area of the tunnel junctions is 
estimated to be ~ 0.01 /im 2 . 

Measurements were done in a dilution refrigerator ca- 
pable to hold temperature from 20 mK to 1.2 K. Mag- 
netic field of 2 Tesla was applied to keep aluminum in 
the normal state. Thermocoax© cables and a copper- 
powder filter next to the mixing chamber of the refriger- 
ator provided necessary filtering against high frequency 
noise penetration. 

We have measured the I—V characteristics of the tran- 
sistors at different temperatures and gate voltages. At 
low voltages our measurements were performed using a 
lock-in amplifier with 6 Hz reference signal frequency in 
the constant voltage mode (a feedback was used to con- 
trol the amplitude of the excitation in order to keep the 
output signal at the same level). 



Theory. At not very low temperatures (or voltages) the 
quantum dynamics of a tunnel junction is well described 
by the quasiclassical Langevin equation O] 
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Jj + £ji cos tpj + £ j2 sin tpj, (1) 



j = L,R; Jj is the current flowing through the junction; 
the phase variable ipj is related to the voltage across the 
junction as fUpj/e = Vj. Discrete electron tunneling is 
responsible for the noise terms in eqs. ([!]) described by 
the stochastic Gaussian variables £jk(t) with correlators 
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The eqs. (Q) are supplemented by the appropriate cur- 
rent balance and Kirghoff equations and can be solved 
perturbatively in the noise terms (see |7],[L0| for details). 
One arrives at the I — V curve for a SET transistor |jj : 
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where Q av (V) = C g V g + HLC B L L +B R Ii Cn V is the average 
charge of the island and i?s = Rl + Rr. The last two 
terms in eq. (g) describe deviations from the Ohmic be- 
havior due to charging effects. In the case Rl = Rr for 
the V g -independent term Iq(V) we find |Q 
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*(x) is the digamma function, a = 1 — ieV/AirT and 
b = a + 2a t E c /Tr 2 T. The expression (|) holds for § 
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7 = 0.577.. is the Euler constant. The last term in eq. 
(3) describes the modulation of the I—V curve by the 
gate voltage. Provided the island charge fluctuations are 
large the amplitude of the modulation is exponentially 
suppressed F(T, V) > 1. The general expression 

for the function F(T, V) [0 is quite complicated and is 
not presented here. In the limit of small voltages and for 
T <C a t Ec this expression becomes simpler, and we get 



F(T, 0) ~ (T/T o y, T Q = VT2^tEc/n 2 



(6) 



As the condition (||) should be simultaneously satisfied 
eq. (||) makes sense only for a t > 1. The constant A 
in eq. (g) has the form HflGj A = f(a,)e- 2at . The 
prefactor f(a t ) was also estimated in [j? 10 



with a limited 



accuracy. More accurate results for f(a t ) at low T can 
be derived by means of other techniques ||. 

Results and Discussion. In what follows we will dis- 
regard a small asymmetry in the parameters of L- and 
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i?-junctions and put Rq = i?s/4. The value of Rs was 
measured from the slope of the I — V curve at high volt- 
ages. The accuracy of these measurements was limited 
by nonlinearities on the I — V curves due to suppression 
of tunnel barriers and heating and is estimated as ~3%. 
The charging energy Ec is usually determined from the 
offset on the I — V curve at large V. In the strong tun- 
neling regime a t > 1 a clear offset can be reached only at 
very high voltages where precise measurements are diffi- 
cult due to other reasons. Therefore, the above method 
gives only a rough estimate of Ec with the accuracy 
within a factor 2. Alternatively one can try to deter- 
mine Ec from the high temperature expansion of a zero- 
bias conductance: G{T) = (l/i?s)(l - E C /3T+ ..). For 
our samples this asymptotics works well only at T > 10 
K and could not be reached in the same cooling cycle. 
Therefore, the parameter Ec was determined from the 
best fit of the low-temperature (< IK) and low-voltage 
(< 700/ieV) parts of the I ~ V curves averaged over the 
gate charge. Fitting of dl/dV for different temperatures 
allows to determine Ec with a sufficiently high accuracy 
and avoid problems discussed above. Also the precise 
value of at can be verified by this method. 

The results are summarized in the following table: 



TABLE I. Parameters of the samples 



Sample 




a t 


E c , K 


10wo, /ieV 


Wwo, mK 


I 


17.4 


1.48 


2.25 


50 


600 


II 


12 


2.15 


1.1 


10 


100 


III 


10.4 


2.48 


1.04 


5 


60 


IV 


6.5 


3.97 


1.16 


0.5 


6 


V 


3.1 


8.32 


~ 0.3 


5xl0 -4 


5xl0" 3 



The last two columns show the lowest voltage and tem- 
perature ~ 10 wo where the theory is still applicable. 

The data for a temperature dependent zero bias con- 
ductance averaged over the gate charge are given in Fig. 
2 for four samples (filled symbols) together with the the- 
oretical dependencies G*y° or = 1/Rs — Iq/V\v=o (solid 
curves). The Coulomb blockade induced suppression of 
the conductance at low T is clearly seen even for a highly 
conducting sample IV. We observe a good agreement be- 
tween theory and experiment except at the lowest tem- 
peratures where the theoretical results become unreli- 
able. In this temperature interval the theoretical curves 
turn out to be closer to the minimum conductance which 
is also shown by open symbols. At low T the system con- 
ductance shows a tendency to saturation. This is com- 
patible with the corresponding conjecture made in Ref. 

. On the other hand, heating effects as an additional 
reason for this saturation cannot be excluded cither. 

The averaged differential conductance dl/dV is shown 
in the bottom panel of the Fig. 3 for all five samples 
at the lowest temperature T s» 20 mK. It is remarkable 
that even for the sample V with a< > 8 a decrease of 
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FIG. 2. The averaged over V g (filled symbols) and the min- 
imum (open symbols) values for the linear conductance of the 
samples I-IV together with a theoretical prediction for the 
averaged conductance (solid curves). 

the differential conductance at small V due to charging 
effects is well pronounced. The values dl/dV measured 
for the sample II for different temperatures are presented 
in the top panel together with a theoretical prediction 
from eq. (^) . A similarly good agreement between theory 
and experiment was also found for the samples III, IV, 
and V. 

The gate modulation of the current was found to be 
of a simple cosine form (||) for all samples except for a 
relatively highly resistive sample I at low T. The results 
for the gate modulated linear conductance of the sample 
II are shown in Fig. 4a. The amplitude of modulation 
increases with decreasing temperature in a qualitative 
agreement with theory. At low T the modulation effect 
is considerable even for the sample IV with at — 4. At 
the lowest temperatures this effect is visible also for the 
sample V. However in this case the amplitude of modu- 
lation was only slightly above the average noise level. 

The data for the temperature dependent amplitude of 
conductance modulation for the samples I-IV are pre- 
sented in Fig. 4b. Solid curves correspond to the best fit 
with a theoretical dependence oc exp(— T 2 /Tq) (cf. eqs. 
([| ||)). Note, that for all samples the best fit value To was 
found to be by a numerical factor ~ 2 3 smaller than 
the value (^). In other words, the temperature suppres- 
sion of the gate modulation is bigger that it is predicted 
theoretically. We speculate that this may be due to an 
additional effect of noise. Another possible reason for 
such a discrepancy is an insufficient accuracy of the the- 
oretical calculations of the gate modulated conductance. 

Our experimental results clearly show that - in ac- 
cordance with earlier theoretical predictions - Coulomb 
blockade is not destroyed even in the strong tunneling 
regime: clear signs of Coulomb suppression of the transis- 
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FIG. 3. The Vg-averaged differential conductance. Top 
panel: the data for the sample II (thinner curves) together 
with theoretical results (thicker curves). For the sake of clar- 
ity the curves at higher T are shifted vertically with a step 
0.02 kfi -1 . Bottom panel: the data for all five samples at 
T ~ 20 mK. 



tor conductance were observed for at as large as 8.3. For 
all at the characteristic energy scale for charging effects 
is set by the renormalized Coulomb gap Eq, however 
such effects remain well pronounced even at T > ££. 
The I — V curves measured for all five SET transistors 
are in a quantitative agreement with the strong tunnel- 
ing theory (?],[l(| , except at very low temperatures where 
this theory is not applicable. Along with the overal sup- 
pression of the conductance its modulation by the gate 
voltage was also observed at sufficiently low T. The mod- 
ulation effect increases with decreasing temperature in a 
qualitative agreement with theoretical predictions. 

Our work was aimed to study SET transistors in the 
nonperturbative strong tunneling regime. In this sense 
our experiment is complementary to that by Joyez et al. 
Q where SET transistors with smaller at were studied. 
Theoretically, a perturbative regime characterized by an 
expansion parameter g = a t /ir 2 <C 1 (see e.g. H,0,§) can 
be easily distinguished from a strong tunneling one with 
exp(— 2a t ) <C 1 (cf. e.g. HQ) except within the interval 
0.5 < a t < 3 where both inequalities are (roughly) satis- 
fied. Combining our results with those of Ref. || as well 
as with the corresponding theoretical predictions flj^fl0| 
we can draw a more definitive conclusion about the va- 
lidity range for both regimes: tunneling can be treated 
perturbatively for at ^ 1 while the nonperturbative tun- 
neling regime sets in for at ^ 2. In the intermediate 
1 ^ oit 2 nonperturbative results can be also applied 
at temperatures not much lower than Eq- 

In summary, we have demonstrated that SET transis- 
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FIG. 4. The gate modulation amplitude G max — G m in of 
a linear conductance as a function of T for the samples I-IV 
(triangles, squares, diamonds and circles respectively). Solid 
lines provide the best fit with a formula A exp( — (T/Tb) 2 ). 
The inset: the linear conductance of the sample II as a func- 
tion of Q g — C g Vg measured at different T. 

tors operate well even if the effective tunneling resistance 
Rq is several times smaller than 6.5 kf2. 
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